Efficiency of higher dimensional Hilbert spaces for the violation of Bell inequalities 
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We have determined numerically the maximum quantum violation of over 100 tight bipartite Bell 
inequalities with two-outcome measurements by each party on systems of up to four dimensional 
Hilbert spaces. We have found several cases, including the ones when each party has only four mea- 
surement choices, when two dimensional systems, i.e., qubits are not sufficient to achieve maximum 
violation. In a significant proportion of those cases when qubits are sufficient, one or both parties 
have to make trivial, degenerate 'measurements' in order to achieve maximum violation. The quan- 
tum state corresponding to the maximum violation in most cases is not the maximally entangled 
one. We also obtain the result, that bipartite quantum correlations can always be reproduced by 
measurements and states which require only real numbers if there is no restriction on the size of the 
local Hilbert spaces. Therefore, in order to achieve maximum quantum violation on any bipartite 
Bell inequality (with any number of settings and outcomes), there is no need to consider complex 
Hilbert spaces. 

PACS numbers: 03.65.Ud, 03.67.-a 
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I. INTRODUCTION 

One of the most astonishing features of quantum- 
mechanics is its nonlocal nature. Separated observers 
sharing an entangled state and performing measurements 
on them may induce nonlocal correlations which violate 
Bell inequalities [H, 0- In contrast, separable states sat- 
isfy all the possible Bell inequalities with any measure- 
ment settings. 

A general setting concerning Bell inequalities is that 
measurements are made on a system, which is decom- 
posed into N subsystems. On each of these subsystems 
one out of m^, i = l,...,iV observables is measured, 
producing ki, z = 1, . . . , iV outcomes each. In almost all 
the cases investigated up to now in order to maximally 
violate them the dimension of the local state spaces of 
the shared entangled state did not have to be larger than 
the number of outcomes of the respective parties. Some 
notable exceptions to it are the bipartite kA = 3 and 
ks — 2 Bell inequalities in Ref. Q , and families of corre- 
lation Bell inequalities with binary outcomes where 
the smallest number of measurement settings was found 
to be niA = 8 and tob = 4. This latter case requires 
states of dimension larger than the number of outcomes 
to obtain maximal violation. 

In the present numerical investigation our aim is two- 
fold. Firstly, wc wish to demonstrate that by including 
marginal probabilities in the Bell inequalities it is further 
possible to reduce the number of measurement settings. 
Then we also show that any bipartite Bell inequality can 
be violated with settings and states in the real Hilbert 
space in the same extent as with settings and states in 
the complex Hilbert space. 
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Actually, we believe that these results arc not only of 
academic interest: On one hand, higher dimensional sys- 
tems have been produced in the laboratory in a number 
of schemes, subjected to Bell-type tests as well. In partic- 
ular in Ref. (3] the experimental violation of a spin-1 Bell 
inequality has been presented using four- photon states, 
while in Refs. Q, Bell-type tests based on the in- 
equality of Collins et al. Q have been performed for or- 
bital angular momentum and energy-time entangled pho- 
tons producing qutrits, respectively. Also, two-photon 
interference experiments have demonstrated time-bin en- 
tanglement up to d = 20 dimensionality On the 
other hand, this investigation can be especially relevant 
in practical applications of quantum information proto- 
cols. For instance, in quantum cryptography the key 
idea is that only local correlations can be created by an 
eavesdropper, thus the only useful correlations must have 
quantum origin. In order to characterize the set of pos- 
sible quantum correlations useful for quantum cryptog- 
raphy applications, it is important to know how effective 
higher dimensional systems are with respect to qubits. 

In particular, in this paper wc considered tight bipar- 
tite two-outcome Bell inequalities corresponding to the 
facets of the convex polytope [3] with up to five set- 
tings 2-89 of Ref. [13], and the 31 cases with up to four 
settings considered by Brunner and Gisin [l^. We note 
that there is some overlap between the two lists. We 
used projective measurements in all cases, since for bi- 
nary outcomes it has been shown [l^ that general POVM 
measurements arc never relevant. The tools used in the 
numerical exploration are gathered in Sec. |TT1 then in 
Sec. mil we give a list of tables presenting the numbers 
corresponding to the maximum quantum violations in 
cases of real and complex qubits (3-dimensional spaces), 
and real qutrits, taking into account degenerate measure- 
ments as well. For all but two inequalities we considered 
such component spaces were sufficient for maximum vi- 
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olation. In one case complex qutrits, and in one case 
real ququarts (4-climensional spaces) were necessary to 
achieve the maximum violation. For both cases the gain 
was marginal, not much larger than numerical uncer- 
tainty. The numbers obtained are discussed in Sec. IIIIl 
and some conclusions are commented in Sec. lIVI Finally, 
in Appendix [X] wc provide a proof on the equivalence of 
real and complex Hilbcrt spaces in reproducing bipartite 
quantum correlations if there is no constraint on the size 
of the component Hilbert spaces. 



II. THE METHOD 

The quantum value of the expression in the Bell in- 
equality is an expectation value of a Hcrmitian operator. 
The maximum expectation value of such an operator is 
its largest eigenvalue. Therefore, to find the maximum 
quantum violation we have to find those measurement op- 
erators for both Alice and Bob whose combination as it 
appears in the inequality gives the largest possible eigen- 
value • This way the parameters to be optimized are 
those of the measurement operators, no parameter of the 
vector enters the problem. The vector can be determined 
as the eigenvector belonging to the maximum eigenvalue. 

As the outcome of each measurement has to be ei- 
ther or 1. the measurement operators to be considered 
arc projectors in the component Hilbert spaces of Al- 
ice and Bob. In case of 2-dimensional Hilbert spaces 
each nondegenerate measurement operator projects to 
a 1-dimensional subspace, which may be defined by a 
unit vector |m) of irrelevant phase as |m)(m|. Such a 
vector can be characterized by 2 parameters, it is con- 
venient to use the two angles on the Bloch-sphere. As 
it turned out to be essential, we also considered trivial, 
degenerate measurement operators as well. Such a mea- 
surement, represented by the zero and the unit operator 
brings always the result and 1 , respectively. Obviously, 
these measurements need not be performed at all, and the 
problem becomes equivalent with a smaller one with less 
measurements. We performed the optimization with all 
combinations of nondegenerate, zero and unit operators. 
For 3-dimensional spaces a nondegenerate measurement 
operator is either a one or a two-dimensional projector. 
A unit vector of irrelevant phase is again sufficient to 
define either a one and a two-dimensional projector as 
|m)(m| and I — \m){m\, respectively. Four real param- 
eters, for example the two polar angles and the phases 
of two components (one component may be chosen real) 
arc needed to characterize such a 3-dimcnsional complex 
vector. Although we have considered only nondegener- 
ate operators, as each of them may be either a one or a 
two dimensional projector, many optimization runs are 
necessary to cover all combinations. In the case of 4- 
dimensional component spaces we confined ourselves to 
2-dimensional projection operators. To make the opti- 
mization of the many parameters involved for all com- 
binations of the dimensions of the operators would have 



taken too much computer time. A 2-dimensional pro- 
jector in a 4-dimensional complex space requires 8 real 
parameters to define. 

We may reduce the number of parameters involved by 
using the fact that both Alice and Bob may choose their 
bases freely. With an appropriate unitary operation we 
may transform one of the operators, say the first one, 
into a diagonal form. This eliminates all parameters of 
that operator. Then we may apply another unitary oper- 
ator that does not affect the matrix of the first operator 
to simplify the matrix of the second operator as much 
as possible. If there exists further transformation that 
leaves the first two matrices unchanged, it may be used 
to reduce the number of parameters of the third opera- 
tor, and so on. Following this recipe, for qubit spaces the 
vector characterizing the first (nondegenerate) operator 
will be one of the basis vectors (no parameter), while the 
one corresponding to the second operator may be trans- 
formed to have both components real (1 parameter). 

In a 3-dimensional Hilbert space the compo- 
nents of a unit vector may be parameterized as 
(cos vJsinT^e*", sin(^sin'!?e*'^, COS!?), with the 3rd compo- 
nent is chosen real (4 parameters). The vector corre- 
sponding to the first measurement operator may be trans- 
formed to (0, 0, 1) (no parameter). This form is invariant 
to a unitary transformation of the ui2 type (operation 
within the subspace spanned by the first two basis vec- 
tors). With such an operation we may eliminate the sec- 
ond component of the second vector, and wc also make its 
first component real, leaving the form (sint92, 0, cos?92) 
(1 parameter). After this we still have the freedom to 
eliminate the phase of the second component of the 3rd 
vector. 

In the case of 4-dimensional Hilbert spaces, the first 
measurement operator may be diagonalized to have the 
form diag{l, 1, 0, 0). Then we may apply a further trans- 
formation of the form U12M34 to simplify the second oper- 
ator. Wc can obviously diagonalize the two 2x2 blocks 
in the upper left and the lower right corners. Then using 
the fact that the matrix corresponds to a 2-dimensional 
projector, it can be shown that the rest of the trans- 
formed matrix must also have a special form, which with 
a further allowed operation may be simplified to the two- 
parameter form of (1 -|-7i)/2, where 1 is the unit matrix, 
and 
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This has been shown in Rcf. 19]. The first two matrices 
leave no further room to simplify the 3rd and any further 
operators, it will take 8 parameters to characterize each 
of them. We have chosen those parameters by using the 
fact that the matrix of the most general two-dimensional 
projector in the 4-dimensional space may be produced 
by applying the most general transformation of the form 
""12^*34 to the two-parameter matrix above. Each of the 
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2-dimcnsional unitary operators U12 and U34 have 4 pa- 
rameters. However, an overall phase is irrelevant, and 
it also turns out that the effect of the transformation to 
the special form will only depend on the difference of two 
phase angles in the operators, which makes it possible to 
eliminate one more parameter, leaving altogether just the 
necessary number of 2 + (2 • 4 — 2) = 8 parameters. 

We determined the maximum violation with both com- 
plex and real Hilbert spaces. A measurement operator in 
the real space needs just half as many real parameters to 
characterize as in a complex space of the same number 
of dimensions. The parameters we used were the same 
as in the complex space with all phase angles taken to 
be zero. For optimization we applied an uphill simplex 
method [3. As such a method climbs to a local max- 
imum, to find the global one we restarted the method 
from random positions many times, at least 10000 times 
for the 4x4 dimensional Hilbert spaces. We still can 
not be sure that we have found all global optima, espe- 
cially for the largest, the 5522 (5 settings of 2-outcome 
measurements for each of the two parties) cases. Nev- 
ertheless, the results calculated with spaces of different 
dimensions are fully consistent with each other. Either 
with complex or real spaces, a higher dimensional cal- 
culation has always given at least as large violation as 
the lower dimensional ones. When we managed to find a 
larger value, some optimization runs still ended up with 
the lower dimensional result. From properties of the opti- 
mum in the higher dimensional case, namely the number 
of terms in the Schmidt decomposition of the eigenvector 
and the relation of the subspace defined by the Schmidt 
decomposition to the measurement operators may reveal 
if it actually corresponds to a lower dimensional case. 
The 4-dimcnsional calculations can and do reproduce all 
lower dimensional results we considered, including the 2- 
dimensional cases with degenerate operators. When the 
Schmidt decomposition shows that the eigenvector oc- 
cupies only 2-dimensional subspaces of Alice and Bob's 
component spaces, and there are measurement operators 
that project to exactly those subspaces, or to their com- 
plementer space, then those measurements for the eigen- 
state do behave like degenerate ones. Actually, we re- 
alized from such analysis that in most cases when we 
found a larger violation with ququarts than with qubits, 
the higher dimensionality was not essential, just degener- 
ate operators had to be considered. In their recent paper 
Brunner and Gisin also concluded that for one of their 
cases they needed degenerate [l^ measurements. The 
4-dimensional calculation reproduces the 3-dimensional 
results too, and may even reveal, which measurement 
operators should be one, and which ones should be two- 
dimensional projectors for maximum violation. 



III. DISCUSSION OF THE RESULTS 

We calculated the maximum violation of the tight bi- 
partite Bell inequalities A2 — Agg listed in Ref. (^1 is 



TABLE I: Maximum quantum violation of Bell inequalities 
calculated with real qubit component spaces, with nondegen- 
erate measurements. Higher dimensional spaces have given 
no larger violation for these cases. Entries when maximum 
violation is achieved by the maximally entangled state are 
marked by stars. 



Case 


Type Qubit (R) 


Case Type 


Qubit (R) 


Urlbrl(/l2 ) 


2222 


0.207107 ^ 


A27 


5522 


r\ A OAT 


T { A \ 


3322 


0.250000 


A 

^28 


5522 


U.D4Uoi4 


r3 

J4322 


4322 


0.436492 ^ 


^30 


5522 




r2 

-'4422 


4422 


0.521371 


A 

A31 


5522 


A IXTOO 1 T 


A 

A5 


4422 


A A COO A 

0.435334 


A 

^35 


5522 


Pi fiOA AAQ 




4422 


0.541241 ^ 


A 

A40 


5522 


A CATOC A 


A C 


4422 


A TO /I no ^ 

0.87o493 


A42 


5522 


A ^ 1 AOi? K 


ATT 


4422 


0.605554 


A43 


5522 


U.DiU / DO 


AII2 


4422 


0.500000 * 


A51 


5522 


0.660781 


J 4422 


4422 


0.436492 * 


At,2 


5522 


0.621861 


J 4422 


4422 


0.461684 


A53 


5522 


0.638610 


7-10 
J 4422 


4422 


0.613946 


A54 


5522 


0.593681 


J 4422 


4422 


0.638354 


At,7 


5522 


0.660344 


t12 
J 4422 


4422 


0.618814 


A5S 


5522 


0.648890 


J4422 


4422 


0.671409 


A72 


5522 


0.696282 


Alo 


5422 


0.415390 


A74 


5522 


0.689069 


A22 


5422 


0.623457 


A77 


5522 


0.665558 


A24 


5522 


0.604799 


A7g 


5522 


0.892702 


A25 


5522 


0.603379 









a trivial 1122 type, which can not be violated). These in- 
equalities are the part involving at most 5 measurement 
settings per party of a huge list of inequalities obtained 
with the method described in Ref. [l^l . We also included 
the 31 known tight inequalities with up to 4 measure- 
ment settings per party considered recently by Brunner 
and Gisin [T^. We adopted the notation used in that 
paper. Out of the 26 inequalities of 4422 type, 20 was 
newly introduced there, while II422 "^^^ presented in [2l| . 
/I422 in [121^5, Aq, Alh and Alh in [H, while ASi 
and AS2 in [S]. The only 2222 one is the Clauser-Horne- 
Shimony-Holt (CHSH) inequality 0]. The Bell inequal- 
ity found in Ref. [2^ is the only 3322 type, and the three 
4322 cases were introduced in Ref. [2l| . The two lists we 
considered have some overlap, we marked those cases in 
our tables. For every inequality in the lists the classical 
value to be violated is zero, except for /4422J where it is 
one. The maximum violations we show in the tables are 
just the maximum eigenvalues we found, except for the 
case 11^221 where it is one less. 

In Table |T] we listed all those cases for which we 
could not find a stronger violation in any of our calcula- 
tions than the maximum violation we achieved with real 
qubits, performing only nondegenerate measurements. In 
all tables we marked with a star the cases when maxi- 
mum violation was achieved with the maximally entan- 
gled state. For most instances this is not so, which has 
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TABLE II: Maximum quantum violation is reached with com- 
plex qubits, no degenerate measurements. 



Case Type 


Qubit (R) 


Qubit (C) 


-^4422 


4422 


0.414214 * 


0.449490 * 


-^4422 


4422 


0.441730 


0.454837 


^8 


5422 


0.555704 * 


0.591650 * 


Ag 


5422 


0.451695 


0.465243 


All 


5422 


0.445211 


0.456108 


Ai2 


5422 


0.452098 


0.487709 


Al5 


5422 


0.447760 


0.449628 


^19 


5422 


0.588932 


0.622630 


A20 


5422 


0.564956 


0.602240 


A23 


5522 


0.528521 


0.546073 


A26 


5522 


0.486495 


0.527555 


A2<3 


5522 


0.456259 


0.492064 


A32 


5522 


0.396861 


0.413553 


A33 


5522 


0.561909 


0.622631 


Ase 


5522 


0.419088 


0.438868 


A37 


5522 


0.456106 


0.486887 


Ass 


5522 


0.428958 


0.469913 


^39 


5522 


0.612269 


0.617203 


A41 


5522 


0.419234 


0.478563 


A47 


5522 


0.402679 


0.460854 


A48 


5522 


0.431439 


0.454841 


A49 


5522 


0.454198 


0.466694 


A50 


5522 


0.500000 * 


0.518290 


A79 


5522 


0.606128 


0.624315 


Asi 


5522 


0.662368 


0.669010 


Ass 


5522 


0.696038 


0.696166 


As5 


5522 


0.610060 


0.641141 


As6 


5522 


0.780438 


0.800443 



also been noted in Ref . [l^ . Table[TTlcontains the inequal- 
ities when we got the maximum violation with measure- 
ments on complex qubits. For the cases in these tables we 
got the same values for maximum violation with complex 
qutrits and complex ququarts than with complex qubits, 
and real qutrits did as well as real qubits. However, with 
real ququarts we could always achieve the same amount 
of violation as with complex qubits. It is generally true 
that if a bipartite Bell inequality with arbitrary outputs 
per party can be violated by a certain amount with pro- 
jective measurements in n-dimensional Hilbcrt spaces, 
than they can be violated by at least as much with projec- 
tive measurements in 2n-dimensional real Hilbert spaces. 
This property is an immediate outcome of an even more 
general statement, which is provided in Appendix [A] It 
is an open question, whether Lemma I A . 1 1 could be some- 
how generalized so that this statement would be true 
for any multipartite Bell inequalities as well. From the 
construction it follows, and we have demonstrated in Ap- 
pendix [XI that the Schmidt-decomposition of the state 



TABLE III; Maximum quantum violation is reached with real 
qubits, with some measurement operators degenerate. 



Case 


Type 


Qubit (R) 


Qubit (C) 


Qubit (R) 






nondeg. 


nondeg. 


deg. op. 




-'4322 


4322 


0.154701 


0.236068 


0.414214 


* 


-^4322 (^4) 


4322 


0.231613 


0.259587 


0.299038 


* 


Ae 


4422 


0.222941 


0.232051 * 


0.299038 


* 


-'4422 


4422 


0.238042 


0.238042 


0.414214 




-'4422 


4422 


0.055979 


0.055979 


0.414214 


* 


rl3 

-'4422 


4422 


0.249466 


0.250000 * 


0.434855 




rl4 

-'4422 


4422 


0.407621 


0.410296 


0.479410 


* 


rl5 

-'4422 


4422 


0.238273 


0.250000 * 


0.434855 




rl6 

-'4422 


4422 


0.240659 


0.240659 


0.414214 


* 


An 


5422 


0.221946 


0.221946 


0.375447 




Ais 


5422 


0.210377 


0.212229 


0.384355 




A34 


5522 


0.461083 


0.513972 


0.535012 




A44 


5522 


0.500000 * 


0.533925 


0.536494 




A55 


5522 


0.451941 


0.486823 


0.621320 


* 


A56 


5522 


0.675426 * 


0.675426 * 


0.689312 


* 


A59 


5522 


0.430220 


0.430220 


0.448826 




Aes 


5522 


0.327627 


0.327627 


0.479410 


* 


Aeg 


5522 


0.330388 


0.330388 


0.609610 




A70 


5522 


0.465198 


0.465198 


0.605223 




An 


5522 


0.418729 


0.418729 


0.449016 




A73 


5522 


0.800326 


0.852797 


0.883138 




A75 


5522 


0.572736 


0.587052 


0.605151 




Aso 


5522 


0.136376 


0.174354 


0.375447 




As2 


5522 


0.314943 


0.314943 


0.454573 




As4 


5522 


0.605340 


0.619437 


0.623457 




Ass 


5522 


0.076842 


0.076842 


0.414214 


* 



TABLE IV: Maximum quantum violation is reached with 
complex qubits, with some measurement operators degener- 
ate. 



Case Type 


Qubit (R) 


Qubit (C) 


Qubit (R) 


Qubit (C) 






nondeg. 


nondeg. 


deg. op. 


deg. op. 


Aie 


5422 


0.416036 


0.416036 


0.446167 


0.457107 * 


A45 


5522 


0.482065 


0.509936 


0.534037 


0.537239 


Aei 


5522 


0.307654 


0.307654 


0.395168 


0.401925 


Ag2 


5522 


0.219048 


0.231812 


0.395168 


0.401925 


Aaa 


5522 


0.345116 


0.360817 


0.452098 


0.487709 



in the 4-dimensional real space has 4 terms, the Schmidt- 
coefficients are pairwise equal, and the ratio of the pairs 
equals to the ratio of the Schmidt-factors from the qubit 
case with the same violation. 

There are surprisingly many inequalities that can be 
violated more, sometimes very significantly more by al- 
lowing measurements to be degenetatc, than by confining 
ourselves only to nontrivial ones. Table Hill and Table HVl 
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TABLE V: Maximum quantum violation is reached with real qutrits. 



Case 


Type 


Qubit (R) 


Qubit (C) 


Qubit (R) 


Qubit (C) Qutrit (R) 






nondeg. 


nondeg. 


deg. op. 




deg. op. 






^4422 (^7) 


4422 


0.197048 


0.197048 


250000 


* 


250000 




2S7S6S 


^4422 


4422 


0.420651 


0.420651 




* 


48431 3 




487768 


TlS 

-'4422 


4422 


0.181236 


0.181236 


543599 




543599 




642967 


rl9 

-14422 


4422 


0.369700 


0.430724 * 


44S587 




443587 




0.497171 


1-20 
-'4422 


4422 


0.305645 


0.305645 


4'^4'^24 




434324 




449669 


Al3 


5422 


0.397412 


0.403098 


0.414214 


* 


0.414214 




419982 




5422 


0.449958 


0.453901 


452465 








464584 


A46 


5522 


0.446602 


0.449849 










0.458105 


Aeo 


5522 


0.252968 


0.252968 


0.375447 




0.375447 




0.390611 


AeA 


5522 


0.375234 


0.375234 


0.375447 




0.375447 




0.390089 


A65 


5522 


0.208545 


0.208545 


0.347759 


* 


0.353146 




0.355021 


Aer 


5522 


0.395696 


0.395696 










0.396289 


Aes 


5522 


0.385731 


0.385731 










0.395718 


A76 


5522 


0.404741 


0.415397 


0.447555 




0.447555 




0.489863 


Asg 


5522 


0.131420 


0.131420 


0.250000 


* 


0.250000 




0.288932 



show the cases when we got the maximum violation with 
real and complex qubits, respectively, taking one or more 
measurements of Alice, or Bob, or of both of them degen- 
erate, i.e., either unity or zero. As wc have already men- 
tioned, the four-dimensional calculations can always re- 
produce these values even by confining ourselves to rank 
2 measurements (2-dimensional projectors) by operators 
that project onto the subspace the eigenvector occupies, 
or onto the orthogonal one. However, when a complex 
qubit result is reproduced with real ququarts, the eigen- 
vector requires the whole component spaces (4 terms in 
Schmidt decomposition), therefore effect of degenerate 
operators can not be simulated with rank 2 operators 
this way. 

Brunner and Gisin [iB l calculated the maximum quan- 
tum violation by applying degenerate measurements only 
for their 14422 inequality. They did that after realizing 
that this inequality can not be violated by the maximally 
entangled state without such measurements. They state 
(l/\/2 — 1/2) as the value of maximum quantum viola- 
tion, which they achieved by taking two measurement 
operators of both parties degenerate. We found twice as 
large maximum violation by taking two measurement op- 
erators of only one party degenerate (see Table . We 
also found that a very small violation may be achieved 
by using only true two-outcome measurement. The vio- 
lating state is far from the maximally entangled state, it 
has Schmidt coefficients of 0.9158 and 0.4016. 

So far we have only shown cases for which maximum 
violation could be achieved in qubit spaces. The exis- 
tence of Bell inequalities for which this is not the case 
has been proved in Refs. 0,0, Particularly, in Ref. Q 
we were able to give concrete examples of correlation Bell 
inequalities (i.e., inequalities without local marginals) 
whose maximal violation is not achieved by qubits. In 



the present list we found numerically quite a few such 
cases, now for Bell expressions with marginals. In all 
such cases except for two, real qutrit spaces were enough 
for maximum violation, sec Tablc|Vl For most of them, in 
two dimensions larger violation can be achieved by allow- 
ing degenerate operators than by not allowing them (no 
entry in the appropriate place, when it is not so). With 
qutrits we can do even better. However, for most entries 
in the list the increase is quite small, no more than a cou- 
ple of percents, sometimes even much less, which means 
these cases may have no practical and experimental rele- 
vance. For a few cases the gain is more than 10%. We find 
the largest increase (about 0.1, or 18%) for 1^22- It is 
interesting to note that there exist Bell inequalities that 
can be violated more with real qutrits than with com- 
plex qubits, and there are also examples for the opposite 
(at least without allowing degenerate measurements for 
qutrits, which we have not tried). For all cases in Ta- 
ble |V] each party has at least 4 measurements, in the 
smallest ones each of them has just 4. We will show in 
a forthcoming publication that for correlation type in- 
equalities to get larger violation with higher-dimensional 
spaces than with qubits, one of the parties must have at 
least 4 measurements, and then the other one must have 
at least 7 measurements. All 4422, 5422 and 6422 cor- 
relation type Bell inequalities can maximally be violated 
by qubits. 

We found one single inequality in the list that we 
could violate more with complex qutrits than with real 
ones or with qubits. The maximum violation of Ai\ 
(5422) with real qubits (no degenerate measurement) is 
0.099090, with complex qubit (no degenerate measure- 
ment) is 0.125000, with real and complex qubit (degen- 
erate measurement allowed) 0.299038 (maximally entan- 
gled state), with real qutrit 0.316523, and with com- 
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plex qutrit 0.317496. The last improvement is absolutely 
marginal, but it does not seem to be due to numerical 
error. 

For Agj (5522) we foimd we need ququarts to get max- 
imum violation, but the improvement was even less con- 
vincing. The best qubit value is 0.756199 (both with real 
and complex qubit), while the maximum we got with 
both real and complex ququarts is 0.756247. From a 
more detailed analysis of the solution we could not see a 
way to reduce it to a lower dimensional space. It turned 
out that this violation could be achieved by taking two 
measurement operators equal. Therefore, we calculated 
the maximum violation with qubits of the 5422 inequality 
we got by uniting these two measurements, and we found 
0.755931, a slightly smaller value than for the original 
inequality. The difference from the ququart value is still 
extremely small, but at least it seems to be more than 
numerical error. 

In our calculations the maximum number of dimen- 
sions for the component spaces were four. Moreover, we 
allowed degenerate measurements only for qubit spaces, 
and confined ourselves to rank 2 measurements in four 
dimensions. For some cases on the list it is possible, that 
without these restrictions one could find a larger maxi- 
mum quantum violation. 



IV. SUMMARY 

Let us briefiy summarize the main results achieved in 
this work. 

We investigated numerically the maximum values on 
tight bipartite two-outcome Bell inequalities in cases 
when the local Hilbert space was restricted to d = 2, 3, 4 
dimensions. We found Bell inequalities with four mea- 
surement settings for each side where qutrits were needed 
to achieve maximal violation, and with five measure- 
ment settings for each side where ququarts were needed 
to achieve maximal violation. We may interpret these 
results via the concept of witnessing the Hilbert space 
dimension 0, @1- The question is that given a joint prob- 
ability distribution of measurement results performed by 
separate parties, is it possible to set a bound on the di- 
mension of the multipartite state space? Thus, dimen- 
sion witnesses are operators 0] able of bounding the di- 
mension of a quantum system. This allows one to test 
experimentally the size of the underlying Hilbert space, 
which otherwise is a rather abstract concept. Therefore, 
by adapting this language, we can say that we found nu- 
merically tight Bell inequalities which act as dimension 
witnesses for qubits and qutrits. 

On the other hand, in analogy to the terminology di- 
mension witnesses one may inquire whether reality wit- 
nesses could be constructed, which would be able to dis- 
tinguish complex Hilbert spaces from real Hilbert spaces. 
Actually, the existence of such kind of a witness has been 
quested by Gisin in Ref. However, according to our 
result presented in Appendix [XI we may safely say that 



reality witness cannot be constructed for the case of two 
parties since by doubling the size of the local complex 
Hilbert space of each party one may reconstruct all the 
joint probabilities with local real Hilbert spaces as well. 
Although, the question is remained open for multipartite 
systems, numerical study supports us to believe that our 
Lemma holds for the most general case as well. 
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APPENDIX A: ON THE EQUIVALENCE OF 
REAL AND COMPLEX HILBERT SPACES IN 
REPRODUCING BIPARTITE QUANTUM 
CORRELATIONS 

Here the following main result is shown: 

Lemma A.l. Joint probabilities between two separated 
observers which has quantum origin can always be repro- 
duced by measurements and states which require only real 
numbers. 

This fact which is interesting by its own, has some 
striking consequences, an immediate one is that the max- 
imum quantum violation of any bipartite Bell inequal- 
ity (with any number of settings and outcomes) can be 
achieved in the real Hilbert space as well. 

To set the scene, we assume that two separated ob- 
servers, Alice and Bob, may perform one of a finite num- 
ber of measurements, and that each measurement has 
a certain number of outcomes. We label outcomes cor- 
responding to different measurements distinctly, so that 
each outcome a and b is uniquely associated to a single 
measurement of Alice and Bob, respectively. Let Sa and 
Sb be n-dimensional complex Hilbert spaces of the two 
parties, respectively, and \V) be any vector in the ten- 
sor product space Sa Sb- Let Pa (A) be projection 
operator associated with outcome a {b) of Sa (Sb)- 

In the light of the above definitions, we say that the 
joint probabilities pab admit a quantum representation 
[2q if there exists a quantum state p on the composite 
Hilbert space, a set of projectors Pq (8) 1 of AHce's and a 
set of projectors 1 (X) Pf, of Bob's system, such that 

Pab = Tr{PaPbp). (Al) 

Note, that since we do not impose any limitation on the 
dimension of the local Hilbert spaces, we may consider 
projection operators instead of the more general POVM 
measurements. The Bell expression consists of a linear 
combination of probabilities (|A1[) . The projectors be- 
longing to different outcomes of a measurement are or- 
thogonal to each other, and they sum up to unity. 

First we prove the following correspondence between 
joint distributions arising from projection measurements 
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in complex n-dimensional local Hilbert spaces and pro- 
jection measurements in real 27i-dimensional local Hilbert 
spaces: 

Lemma A. 2. There exist projection operators and 
P^ of the 2n- dimensional real spaces S'j^ and S'g, respec- 
tively, and \V') € 5^ ® S'g such that the corresponding 
expectation values are equal, i.e., 



{V\Pa®Pb\V) = 



(A2) 



where the state \V) and operators Pa, Pb cire defined 
above, and \ V'), P^ and P^ depend only on \V) , Pa and 
Pb, respectively. 

Proof. Let us use a matrix representation. Let us choose 
orthonormal bases in each component space, and let the 
basis in the product space be the basis consisting of the 
products of the basis vectors of the component spaces. 
Hence, we can write. 



\V)=J2v.,\vf)\v' 



and 



\Pb\v 



(A3) 



(A4) 
(A5) 



where the basis vectors {Iw/')}"^! and span 
respectively Alice and Bob's local state spaces. This 
way the vectors of the product space will be represented 
by matrices of two indices. Then the expectation value 
above can be expressed as 

V*AkBjiVki=TriAVB^V^) 

where A, B and V are the matrix representations of Pa 
and Pb and | V) , respectively. The value of the expression 
is a real number, as it gives the expectation value of a 
Hermitian operator in the product space. 

Let us consider the following mapping Let 
us replace each component Vi ^ vf' + iv( of the n- 
dimensional complex vector with the two-element real 
block of {vf-,vl), and each component Aij = Af^ -|- iAl^ 
of a two-index matrix with the 2x2 block of 




One can prove that the image of the product of cither a 
matrix and a vector, or two matrices will be equal to the 
corresponding product of the images. For n ~ 1 this is 
easy to show. For n > 1 the multiplication in the 2n- 
dimensional space may be done block-by-block, yielding 
the correct result. The mapping also conserves the linear 
combinations of both vectors and matrices. When trans- 
posing matrices one has to be careful. The image of the 
transpose of a matrix will be the transpose of the image 



of the complex conjugate of the matrix. The complex 
conjugation is needed to get the 2x2 blocks right (they 
are not transposed in the image of the transpose). Her- 
mitic conjugation is preserved by the mapping. It is also 
easy to see that the trace operation on the image will 
give a real number, which is twice the real part of the 
value calculated for the original complex matrix (in each 
block the real part of the diagonal matrix element will oc- 
cur twice, while the imaginary part will be off-diagonal). 
Given these rules in hand it is clear that the image of 
a projector is also a projector, the images of orthogonal 
projectors are orthogonal projectors, and if matrices sum 
up to unity, their images will do so, too. Therefore, the 
images of a set of measurement operators will satisfy the 
properties required. 



Let \V') be the vector in S', 



S'n whose matrix V' 



is constructed with the above rule for 2-index matrices 
from the matrix V of \V), and then multiplied by l/\/2 
to get it properly normalized. We note that the mapping 
rule to be applied in the product space is not the same 
as the one applied in the component spaces. That rule 
would actually give just 2n^ components instead of the 
(27i)^ ones. Let there be the matrix of and P^, i. e. 
A' and B' the image of A and B*, respectively. Then 
A'V'B''^V"' will be the image of {1/2)AVB^V'< , the fac- 
tor of 1/2 is occurring due to the 1/V2 normalization 
factor in the construction of V from V. As the trace of 
AVB^V'^ , which is the expectation value in the complex 
space, is real, its value is one half of the trace of its im- 
age, i.e., it is equal to the trace of A'V' B'^V''^ , which is 
the expectation value in the real space. 

□ 

Note that for an arbitrary mixed state p ~ '^2 \i\Vi) (Vi] 
the expectation value Tr(PaPbp) is the convex sum of the 
expectations (|A2[) with coefficients A^, which entails the 
main result Lemma lA. II we wanted to show. 

Aside from its conceptual interest, we mention two in- 
teresting situations where this fact may prove to be use- 
ful beyond justifying our numerical experience that real 
ququarts could yield at least the same amount of viola- 
tion as complex qubits. 

On one hand, in the inequality presented by 
Bechmann-Pasquinucci and Gisin in Ref. [3| having three 
and two measurement outcomes per Alice and Bob, 
respectively, the maximum quantum violation can be 
achieved with projective measurements sharing a max- 
imally entangled state of dimension 3. However, numer- 
ical evidence suggests that using measurement settings 
which require real numbers, the optimum quantum viola- 
tion could not be reached. It arisen as a natural question 
[2^ 1 whether a higher value could be achieved by using 
only real numbers but allowing to occupy larger Hilbert 
spaces. Our result gives the answer in negative for this 
question regarding this particular Bell inequality and also 
prove conclusively that all bipartite Bell inequalities can 
be maximally violated by quantum states and measure- 
ment settings which need in an appropriate basis only 
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real numbers. This latter problem for the general mul- 
tipartite case was posed by Gisin (see also problem 32, 
fundamental questions number 11 in Ref. [28|). 

On the other hand, in Ref. [2^ a hierarchy of con- 
ditions has been formulated through a semidefinite pro- 
gram [29| . This approach can be used for instance to 
obtain upper bounds on the quantum violation of arbi- 
trary Bell inequalities. In this case, however the matrix T 
in question, which should satisfy the positive semidefinite 
constraint is in general Hermitian. Our results, however, 
entails that this matrix needs to be in fact real valued, 
i.e., must be a symmetric matrix. This stronger condition 
thus may provide us with a tighter upper bound on any 
bipartite Bell inequality, than the one which originally 
required the weaker Hermitian condition. 

Now let us illustrate with a simple example, consisting 
of a qubit at each party, the method how to obtain the 



projection operators and the respective states from the 
original complex valued ones. In this case the state of two 
qubits can be written in an appropriate basis as \V) ~ 
a\vi)\vi) + f3\v2)\v2), where the a and (3 Schmidt coef- 
ficients are non- negative numbers, their square adding 
up to 1. Thus the matrix V in Eq. (jA3|) takes the fol- 
lowing simple form, diag{a, (3) whereas a non-degenerate 
projector on the state space of Alice and Bob can be 
written as P^, = (1 ± va)/2, ly £ a,b. Applying the 
mapping rule, discussed above, we obtain the following 
real valued 4x4 matrices, V = {l/^/2)diag{a,a, P, P), 
implying the entangled state (with nonzero a and /3) in 
the 4-dimensional state space, \V) = (a|00) -|- a|ll) -I- 
/3|22) -I- P\33))/y/2 and the corresponding projection op- 
erators = (1 ± i7(T')/2, V G a,b, where a'^ = ^ 1, 
a' — —(Ty ® ay and a'^ = a^® 1- 
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